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The exact renormalization group method is applied to many-fermion systems with 
short-range attractive forces. The strength of the attractive fermion-fermion interaction 
is determined from the vacuum scattering length. A set of approximate flow equations 
is derived including fermionic bosonic fluctations. The numerical solutions show a phase 
transition to a gapped phase. The inclusion of bosonic fluctuations is found to be signif- 
icant only in the small-gap regime. 
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Attractive forces between fermions play a crucial role in many areas of many- 
body physics. In a system of fermions, any attraction can cause such particles to 
form correlated "Cooper pairs" so that a phase transition occurs. The order param- 
eter for this is the energy gap, which appears in the fermion spectrum. Within such 
systems we can identify a weak-coupling regime, with no two-body bound states, 
which corresponds to Bardeen-Cooper-Schrieffer (BCS) superconductivity and a 
strong-coupling regime, with a deeply bound two-body state, which manifests itself 
in a form of Bose-Einstein Condensation (BEC). 

A new and promising way to treat many-fermion systems with these different 
regimes is provided by the Exact Renormalization Group (ergP. This has been 
successfully applied to a variety of other systems in particle and condensed-matter 
physics ■ . The goal of the ERG approach is to construct the Legendre transform of 
the effective action, T, which generates the 1PI Green's functions. One introduces 
an artificial renormalisation group flow for T, depending on a momentum scale k, 
and defines a running effective action by integrating over components of the fields 
with q < k. The RG trajectory then interpolates between the classical action of the 
underlying theory (at large k), and the full effective action (at k — oj^. 

Here we study a system of fermions interacting through an attractive two-body 
contact potential. We take as our starting point an EFT that describes the s-wave 
interaction of two fermions with scattering length arj. Since we are interested in the 
appearance of a gap in the fermion spectrum, we introduce a boson field whose vac- 
uum expectation value (VEV) describes that gap and so acts as the corresponding 
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order parameter. We take the following Ansatz for T 



(1) 



Here M is the mass of the fermions in vacuum and m is taken to be 2M. 

We expand the potential U about its minimum, </>t = (p — A where A is the 
fermion energy gap. To quartic order this has the form 

U(<p, 0t) =UQ + Ul (0t0 _ A 2 ) + | U2 (^»t^ _ A 2 ) 2 . (2) 

In the RG evolution, we start at high k from an almost free bosonic action where 
the VEV of c/> is zero. When k is lowered, we expect that at some scale the system 
undergoes a transition to a phase with spontaneously broken U(l) symmetry and a 
nonzero fermion gap A. 

The evolution equation for T in the ERG has a straightforward one-loop 
structure^ and in the case of constant fx it has the form 



d k r = -*-Tr 



(r^, - RbT 1 d k R B \ + \ Tr [(r^ - Rf)- 1 d k R F \ • (3) 



(2) 

Here ^ F f(bb) is the matrix of the second derivatives of the effective action with 
respect to the fermion (boson) fields and Rb(f) is a matrix containing the corre- 
sponding boson (fermion) regulators. The usual mean-field results can be recovered 
if we include fermion loops only (the second term of this equation). By evaluating 
the loop integrals in we get a set of coupled equations for the coefficents u n and 

at constant [i. Details can be found in our paper** 1 . The initial conditions are 
obtained by assuming that in vacuum our theory reproduces the scattering length 
do, and that the necessary subtractions are identical in matter and in vacuum. 

Although our results can be applied to many systems, for definiteness we 
concentrate on parameter values relevant to neutron matter: M = 4.76 fm _1 , 
Pf = 1-37 fm _1 , and large two-body scattering lengths (|ao| > 1 fin). We have 
checked the dependence of our results on the starting scale K and find that this 
is undetectable as long as K > 5 fm" 1 . Some typical solutions for the evolution 
equations are given in Fig.^ for the case of infinite oo- We compare two different 
approximation schemes, one where we keep fermion loops only, and one where we 
include boson loops as well and we allow to run. For large values of k the system 
remains in the symmetric phase. At k ~ 1.2 fm _1 , the first derivative of U at = 
vanishes and the system enters the broken phase. 

We find that the contributions of boson loops are small. Indeed in the symmetric 
phase they have essentially no effect. Below the transition they do become visible, 
particularly in U2- However their effects on the gap are even smaller, at most ~ 1% 
for 1/\pfclq\ < 1. In this region, the main effect of the bosons is a small enhancement 
of the gap related to the reduction in ui. 
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Fig. 1. Numerical solutions to the evolution equations for the relevant parameters for infinite 
arj and pp = 1.37 fm, starting from K = 16 fm — 1 . The black line shows the full evolution with 
both fermion and boson loops; the grey line the result with fermion loops only. All quantities are 
expressed in appropriate powers of fm — 1 . 

We have also examined the behaviour of the gapover a wider range of 1/(j>f»o)- 
The overall picture is the same as in Marani et al!^ for fermion loops only, with a 
crossover from BCS pairing (with // ~ ep) for l/(ppao) < to BEC (with \i < 0) 
for 1/ippao) > 0. Deviations from mean-field behaviour are present in the BEC 
region and become increasingly noticable for weaker couplings or lower densities. 

In the case of neutron matter, we find gaps of the order of 30 MeV. There is 
a simple explanation for such unrealistically large valued^. For weak couplings, the 
gap is proportional to exp (— (7r/2) cot8{pp )). For nucleon-nucleon scattering, cot 5 
increases quickly with momentum, resulting in a reduction of the gap. We therefore 
expect that inclusion of the effective range should capture this physics. 

There a number of improvements which should be made to our approach. Adding 
an effective range is clearly an important one. Notes that this will also require 
inclusion of the three-body forces to respect the reparametrization invarianc M. The 
running of the boson kinetic mass (Z m ), the fermion renormalization constants and 
the "Yukawa" coupling are needed for a full treatment of the action (JJJ. Beyond 
that we would also like to treat explicitly particle-hole channels (RPA phonons). 
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